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ABSTRACT 



Oh 



I I Context. In contrast to the more extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. 
.^^ Aims. We study the molecular gas in the Galactic high-mass star-forming region NGC 63341, which contains diffuse, quiescent components that 
■ are inconspicuous in widely used molecular tracers such as CO. 

Methods. We present HerschelfHIFl observations of methylidyne (CH) toward NGC 63341 observed as part of the "Chemical HErschel Surveys 
of Star forming regions" (CHESS) key program. HIFI resolves each of the six hyperfine components of the lowest rotational transition (7= \-\) 
of CH, observed in both emission and absorption. 
*Y^ Results. The CH emission features appear close to the systemic velocity of NGC 63341, while its measured FWHM linewidth of 3kms"' is 
Q smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is ~7xl0"", two to three orders of 
?H magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC 63341, we do not detect any 
^ outflow signatures in CH. At least two redshifted components of cold absorbing material must be present at -3.0 and +6.5 km s"' to explain the 
d absorption signatures. We derive a CH column density (A'ch) of 7x 10'^ and 3x10''' cm"- for these two absorbing clouds. We find evidence of two 
additional absorbing clouds at +8.0 and O.Okms"', both with A'ch s;2x10" cm"^. Turbulent linewidths for the four absorption components vary 
between 1.5 and 5.0km s"' in FWHM. We constrain the physical properties and locations of the clouds by matching our CH absorbers with the 
absorption signatures seen in other molecular tracers. 
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^ Conclusions. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas 

O" the CH molecule traces more quiescent material. The four CH absorbing clouds have column densities and turbulent properties that are consistent 

^^ with those of diffuse clouds: two are located in the direct surroundings of NGC 6334, and two are unrelated foreground clouds. Local density and 

^'i dynamical effects influence the chemical composition of the physical components of NGC 6334, which causes some components to be seen in CH 

^~~^ but not in other tracers, and vice versa. 
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r| 1. Introduction ( [Rydbeck et al.|1973[ [Turner & Zuckerman|1974| . The accurate 

. J^ determination of a CH abundance from its radio lines is com- 

^ Most studies of star-forming regions focus on dense gas in pUcated by population inversion effects, while its pure rotational 

/S the cores, while much less is known about surrounding diffuse transitions lie in the far-infrared and are unobservable from the 

^ clouds. Methylidyne (CH) can be used as a tracer of low density ^^^^^ ^j^ij^ tj^^ rotational spectrum of CH has been studied 

clouds, since its column density is known to coirelate with the ^-^^^ ^^e 1970s (see David son et al. | |200T) , the first astronom- 

rf ^, T?Sx^° T/^"''l'^. o^'^'l ,' ^y^o'^'!^"^':''' ical detection of a rotational transition of CH (at 149 /zm) was 

&Bk^kl986}|niF&LS3is 2002; Sheffer etal. 2008 ; Chastain ^.^ .^^^ by .Stacey et al.| ( fT987) in Sgr B2. With the Heter odyne 



et al.|2010f. In conttast, most carbon is locked Tn CO in higher instrument for the Far-llifrEid (HIFI, |De Graauw et alllSoTO^ 



density environments (^(Ha) > 10^ crn'^), resulting in a lower ^^^^^^ ^j^^ i/er.c/ze/ Space Observatory dPilbratt et al.|2()10l, it 



fractional abundance of CH (e.g.,|Polehampton et al.|2005|. -^ ^^^ ^^^^-^^^ ^^ ^^^^^rve the ground-state rotational transition 



Interstella r CH was first observed in its electronic transitions ^f ^H at 533 and 537 GHz, free from atmospheiic interference 

in the visual ( lAdams|1941K and later the A-doublet at 3.3 GHz ^^^ ^^^1^ unprecedented spectral resolution. 

* Herschel is an ESA space observatory with science instruments I" this Letter, we study spectrally resolved signatures of 

provided by European-led Principal Investigator consortia and with im- the CH ^Ua/i J -3/2-1/2 transition, with an upper level energy 

portant participation from NASA. of 26 K, toward NCG 63341, which exhibits emission from the 
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hot core itself, and absorption due to cold foreground clouds. 
NGC 63341 is a high-mass star-forming core at a distance of 
1 .7 kpc belonging to the luminous NGC 6334 complex of molec- 
ular clouds. The "I" source is composed of several dense sub- 
cores ( [Hunter et al.||2006l l with a total m ass of ~2QQM n and 



a bolometric luminosity of several 10^ Lq (Sandell 2000|. It is 



rich in molecular lines (e.g., [Schilke et al. ''2006). Several fore 
ground clouds have been identified along the line of sight toward 
NGC 6334 by [Brooks & Whiteoak[ ( [25oT] ) and [Beuther et al 
( [2005| l. We use Herschel/HIFI data to estimate the CH abun- 
dance in the NGC 63341 hot core and to constrain the physical 
properties of the CH absorbing foreground clouds. 

2. Observations and data reduction 

We use data from a spectral scan from the "Chemical HErschel 
Surveys of Star forming regions" key program (CHESS; 
Ceccarelli et al.|2()T0 1, observed by HIFI on board the Herschel 
Space Observatory. The band la scan was performed on 
February 28, 2010 using the dual beam switch (DBS) mode with 
continuum stability optimization. The half-power beamwidth 
at the relevant frequencies is about 42" and was centered at 
a'/20oo=17''20"53':32, Jyzooo =-35 °46' 5 8.5 Vwhichjs between 
the SMAl and SMA2 cores in Fig. 1 of [Hunter et"aL] ( [2006| . 
The reference positions for the DBS were ~3' east and west of 
the target. HIFI spectral scans use the wideband spectrometer 
(WBS) with a frequency resolution of 1.1 MHz, corresponding 
to 0.6 km s' at the frequencies of the two groups of CH lines. 
The full band la (480-560GHz) scan took 12 950 seconds, in- 
cluding overheads. The resulting noise level is 12 mK in r^b 
units, measured in 0.5 MHz channels, in the 'single sideband' 
spectrum (see below), averaged between H- and V-polarization. 
The double sideband system temperature varied between 60 and 
llOK. ^^-^ 

The raw data is pipeUned in HIPE 2.6 ( [Ott|2010| l. We note 
that the absolute velocity scale is uncertain because of improper 
calculation of the spacecraft motion, but this has been con- 
strained to 0.3 km s"' by subsequent reprocessing in a more re- 
cent pipeline version. The resulting level 2 product was pro- 
cessed in CLASSn After mixing the sky signal with a local os- 
cillator (LO) signal, a heterodyne detector such as HIFI yields 
a double sideband spectrum. In spectral scans, each frequency 
channel is observed using several LO settings (observations dis- 
cussed here use a redundancy of 8), which makes it possible 
to mathematically disentangle the signal from the two separate 



sidebands (|Comito & Schilke 



2002 1. This 'deconvolution' step 



is executed in CLASS, while allowing channel-to-channel varia- 
tions in receiver gain. Resulting gain variations are 10% at max- 
imum. Spectra presented in this Letter are always deconvolved 
single sideband spectra. Intensities are corrected for a main beam 
efliciency of 70.8% expected for band la of HIFI. 

3. Line profile analysis 

The single sideband HIFI WBS spectra near 536.76 and 
532.73 GHz (Fig.fl]) show continuum emission with Tm<o-\.\9 K 
and several spectral lines. The two spectral segments, each about 
0.2 GHz wide, cover the CH ^Hs/a 7=3/2-1/2 parity + to " transi- 
tions around 532.73 GHz and the parity to + transitions around 
536.76 GHz (Davidson et al.,,2001) . Most lines are due to CH, 

' Version of April 2010, with some improvements. CLASS 
is part of the GILDAS software developed at IRAM: 
http : //www . iram.fr/IRAMFR/GILDAS. 
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Fig. 1. Observed spectra (black) of the two groups of three 
hyperfine components of the CH ^113/2 7=3/2-1/2 transition, 
represented in T^\, units along the vertical axis (see Sect. l2]l. 
The horizontal axis gives frequency (top) and velocity with re- 
spect to the rest frequency for the strongest line in each triplet 
(bottom). The magenta line represents model I with one emis- 
sion component and two absorption components and the blue 
line represents model II with two additional absorption compo- 
nents. Physical parameters of the model components are given 
in Table [T] 'Comb' labels indicate the loci of the CH hyperfine 
lines shifted with respect to the rest frame frequency. 



Table 1. Model component parameters. 
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but two emission features in the 532.73 GHz spectrum have con- 
tributions from the CH3OH 7;f=152-15i and the '^CHsOH 93- 
92 transitions. 
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3.1. Model setup 

From inspection of the observed spectrum in Fig. [T] we note 
that the two sets of three emission peaks are spaced by exactly 
the hyperfine splitting. Since we see more than three absorption 
features near each emission triplet, we conclude that more than 
one absorption component is necessary to explain the observed 
spectra. We therefore generate a model CH spectrum starting 
from the minimum of one emitting and two absorbing compo- 
nents ('model r, see Table fl]) using the CASSIS softwara^and 
the CDMS spectroscopic database ( Miiller et al. 2005 1. Given an 
excitation temperature and a column density for a set of emission 
and absorption lines, the 'LTE' module in CASSIS calculates a 
frequency-dependent optical depth, which is used to construct a 
beam-averaged, frequency-dependent line brightness. 

We adopt a source size of 10" for the emission component, 
taken to encompass the four SMA cores ( Hunter et al.| 2006i. 
With a beam size of 42" for HIFI band la, this results in a beam- 
filling factor of 6%. Since we lack collisional rate coefficients for 
CH, we take the dust temperature of 1 00 K from |Sandell| ( |2000| l 
and assume this to be equal to the excitation temperature (Tex) 
of the emitting CH gas in NGC 63341. We note that, while we 
assume a single value for Tex, the nominal dust temperatures 
determined from interferometric observations by [Hunter et al. 
( |200 6 ) vary between 33 and 100 K from core to core, with lim- 
its of 20 and 300 K. The continuum level for the model is fixed 
at the observed level of 1.19 K. For the absorbing foreground 
components, rex=5.5K is based on H2O absorption signatures 
discussed in Emprechtinger et al. (20I0]l. 



3.2. Column density estimates 

Since all features in the 536.76 GHz spectrum (Fig. [T] top 
panel) can be ascribed to CH, we use this as a starting 
point to constrain the remaining model parameters for the 
three components: line of sight velocity (Visi), column density 
(A^ch), and full-width at half-maximum of the line (FWHM). 
The best-fit set of parameters presented in Table [T] (model I) 
are: emission at yisr=-8.3kms"' with FWHM=3.0kms"' and 
A^CH=2.0xlO'^ cm"^ and two absorbing components at -3.0 and 
H-6.5kms ', with distinctly different line widths of 5.0 and 
1.5 km s', respectively. After fine-tuning the parameters of the 
three components, some discrepancies between model I and the 
observed spectra remain. Extra absorption appears to be required 
near 536.782, 536.748 and 536.732 GHz (see Fig.[T] top panel). 
The frequency spacing between these discrepant regions is such 
that they cannot be explained by a single additional component. 

The discrepancies can be resolved by adding two absorbing 
clouds, both with rex=5.5K and A^ch=2x10'^: component C 
with yisr=-i-8.0kms"' and FWHM=3.5kms"', and component 
D with Vki-^O.Okms-i and FWHM=3.0kms-i (Fig. [!}. Upon 
introducing components C and D, the column densities of the 
absorbing components A and B need to be slightly decreased to 
6.6x10'^ and 3.0x10'^ cm"^ respectively ('model IT, Table[T}. 

Without specific tailoring of the model to match the 
532.73 GHz spectrum, it explains this second set of CH signa- 
tures remarkably well, with the exception of the contamination 



ept 




by CH3OH and '^'CHsOH (see Fig. [Ij bottom panel). Although 
the addition of components C and D moves the model closer 
to the spectrum of the 536.76 GHz triplet, the correspondence 
to the 532.73 GHz spectrum still has some notable discrepan- 



cies: the absorption feature near 532.78 GHz is not deep enough, 
and the absorption signature near 532.72 GHz appears to be 
shifted by ~2kms"'. The latter mismatch could be alleviated 
by adding a methanol emission line. Additional contamination 
by other species is not ruled out at the level of ~0.05 K line in- 
tensity, as illustrated by emission features around both 532.670 
and 532.754 GHz, and possibly the blue wings of the methanol 
lines. 



3.3. Uncertainties 

Variations in column density for an individual model component 
have a roughly linear effect on the equivalent width of the corre- 
sponding spectral feature, i.e., doubling A^ch (dividing by two) 
leads to an emission or absorption line that is twice as strong 
(weak). Moreover, the column densities that we derive for each 
component depend directly on the assumed values of Tex. We 
briefly explore the effect of changing Tgx on the line intensi- 
ties, which we compensate by adjusting the column density. For 
the emitting gas, raising Tex from 100 to 150 K requires a com- 
pensating column density increase of ~35%, while lowering Tex 
from 100 to 50 (to 30) K requires a column density decrease of 
only 20% (15%). Because the CH transition under consideration 
becomes optically thick at low Tex, the dependence of column 
density on Tex weakens and eventually reverses as Tex is low- 
ered to ~50K and beyond. To illustrate, if Tex is varied from 150 
to 100 to 50 to 30 K, the opacity for the main emission line in- 
creases from 0.05 to 0.07 to 0.19 to 0.39, respectively. The line 
opacity therefore reacts more strongly than linearly to a change 
in Tex, resulting in counterintuitive behavior of line strength as 
a function of Tex. For the cold absorbing gas, an increase in Tex 
from 5.5 to 20 K roughly doubles the column density, and chang- 
ing Tex from 5.5 to 2.7 K results in a column density decrease of 
<10%. This illustrates a much weaker dependence of column 
density on Tex than for the emitting gas. Optical depths reach 
maxima of 0.34, 0.51, 0.15, and 0.18 at the respective line cen- 
ters for absorbing components A, B, C, and D. 

The uncertainty in our derived column density for the emis- 
sion component is dominated by the uncertainty in the Tex es- 
timates (< 2, see Sect. |3.1| l. Following the above discussion, 
quoted column density values can be trusted taking into account 
an error margin of ~30% for the warm emitting gas. On the other 
hand, the uncertainty in column densities for the cold absorbing 
gas components is dominated by measurement errors in the con- 
tinuum and absorption line depth, each estimated at 6%, result- 
ing in an uncertainty of 8% for the column densities. 



4. Discussion 

4.1. CH in the hot core 

The velocity of the emission component at -8.3 km s"' is con- 



sistent with the H2O emission detected by Emprechtinger et al. 
( ,2010) at -8.1 kms"' and close to the velocity of * -7kms"' 
of mo lecular lines measured from the ground ( |Nummelin et al 



19981. The width of the CH emission features of 3.0 km s ^ is 



smaller tha n the 4-5 km s ' for NH3 (ATCA at 24 GHz; |Beuflier 
et al.|2005| , 6.5 km s' for SiO, and 6kms-' for HC3N (IRAM 



2 CASSIS has been developed by CESR-UPS/CNRS 
(|http://cassis.cesr.fr^. 



30m at 87, 130 and 155 GHz; [Bachiller & Cernicharo|[T990l l. 
This difference suggests that the CH emission originates in dif- 
ferent, less turbulent and/or less evolved, cores than the heav- 
ier N- and O-bearing molecules. Lis et al. ( |2010 i find a sim- 
ilarly narrow linewidth of 3.3kms~' from HIFI obse rvations 
of H^^Cl 1-0 emission from the NGC 63341 hot core. ISandelll 
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( 2000| l derives a source size of (10x8)" and a dust temperature 
of 100 K from sub-millimeter dust continuum measurements. 
Combining this with a 350 //m flux from SHARC-11, dust grain 
opacities from |Ossenkopf & Henning p994| l and the distance of 
1.7kpc, we arrive at A^Hj^SxIO^"* cm"^. From this A^h, and our 
A^CH=2.0xlO''^cm-2, we find an abundance for CH of 7x10 " 
relative to H2 in the hot core. This is a factor ~7-30 lower than 
the CH abundance for the dense envelope of Sgr B2, and a fac- 
tor ~300 lower than observed CH abundances in diffuse clouds 
( Polehampton et al.|2005 Table 4, and references therein). This 
decrease in CH abundance is similar to that predicted by chemi- 
cal models (Sect.[T]l for H2 densities of >10"^ cm" 



,-3 



4.2. CH in absorbing clouds 

The range of turbulent linewidths for our four absorbing clouds 
(Table [1]) and A^ch values of several lO'^ cm"^ agree with mod- 
els of 'turbulence dissipation regions' (Godard et al. 2009|l 
and with optical observations of diffiise clouds ( [Liszt & Lucas| 
|2002 ). Using the range of A^ch values for ab sorption clouds A-D 
(Table[T|l and an A^ch/A^h, ratio of 3.5x10"^ ( [Sheffer et al.|2008| l, 
we find A^H,=(0.6-2)xl0'2' cm^^. Dividing A^Hj of 10^^ cm"^ by 
a number density of 10^-10^ cm"^ applicable to diffuse gas, we 
derive a rough linear size estimate of 0.3-3 pc for the absorbing 
clouds. We also convert A^h^ =(0.6-2)x10^' cm"^ into a visual 
extinction, Ay ~ 0.3-1, based on A^hMv=1-87x10^' cm^^mag"' 
(Bohlin et al. 1978[ l. Hence, we conclude that the CH absorp- 
tion components are diffuse interstellar clouds. We note that, 
although we discuss several absorbing components as separate 
entities in velocity space, a diffuse interstellar cloud of several 
parsecs in size can enjoy a velocity gradient of a few kms"' 
(cf . FWHM in Table [Til and therefore various components may 
originate in a single molecular cloud. 

Other studies also find absorbing gas along the line of sight 
toward NGC 63341. We summarize those that potentially cor- 
respond to our CH clouds. [Brooks & WTiiteoak ( 2001[ l note 
two clouds absorbing in OH (near 1.6 GHz) toward several 
NGC 6334 cores: one cloud with velocities ranging between -15 
and -H2kms"', which they assume to be associated with the 
NGC 6334 complex, and another at -1-6 km s' which they as- 
cribe to a more extended foreground cloud. In addition, the HlFl 
study of H2O lines toward NGC 63341 by jEmprechtinger et al. 



(12010) reports absorbing foreground clouds at -6.3, -0.1 and 
H-6.1kms"'. 

In terms of V^r, our absorption components A and D 
(Tablefljl fit in the widespread [-15, +2] kms ' OH absorption 
clouds, which places them physically close to the NGC 6334 
molecular cloud. In addition, cloud D coincides with one of the 
three H2O absorbing clouds. Hence, although A^ch for cloud A 
is three times higher than for cloud D, we conclude that cloud 
D must be of somew hat higher number densi ty to allow for H2O 
formation. Moreover, [Ossenkopf et al. [ ( [20 1 0[ ) find H20^ absorp- 
tion toward NGC 63341. There is notable uncertainty regarding 
the rest frequencies of the H20^ ground state rotational tran- 
sitions. We assume a value of 1115.209 GHz for F=5/2-3/2, 
which is close to predictions by Miirtz et al. ( 19981), to derive 
Visi~-1 kms ' for the H20^ absorption. Both the above H20^ 
signature and the detection of H^Cl^ abso rption at -1.7kms ' 
with a width of 12 km s"' ( Lis et al.|2010 i are consistent with a 
blend of our CH components A and D in velocity space. 

Our component B at -H6. 5kms~^ is consiste nt with both the 



We believe that this cloud has a relatively high particle density, 
but is unrelated to the NGC 6334 complex. Cloud C is seen in 
CH but has no counterparts in H2O or OH. The column density 
of this cloud is not particularly low with respect to clouds A, B, 
and D (cf. Table [T]), but it might simply be of lower number den- 
sity. In conclusion, from the correspondence with counterparts 
in other molecules, we find that CH clouds B and C are unre- 
lated foreground clouds, while clouds A and D are located in the 
direct surroundings of NGC 63341. 

4.3. Dense clouds and outflows 



In addition to the absorbing clouds discussed in Sect. 4.2 the 
literature provides a collection of absorption signatures with- 
out CH counterparts. First, we find no evidence of absorbing 
CH gas that matches the -6.3 kms ' component seen in H2O 
( Emprechtinger et al. 2010). The low CH abundance in this cloud 
- close to the NGC 63341 hot core in velocity space - may be 
due to its physical proximity to the luminous NGC 6334 region, 
where higher gas densities facilitate active chemistry, leading 
to depletion of CH. In general, overall column density and vi- 
sual extinction are traced more accurately by CH than by H2O 
and CO, which are very sensitive to inhomogeneity (e.g., Spaans 
n^96| ). 

Second, [Beuflier et aT] (p005 ) detect CH3OH and NH3 ab- 
sorption at 25 GHz toward NGC 6334-I-SMA3 (the compact H 11 
region) at -9 kms"', about 2 kms"' blueward of the emission 
component. Moreover, outflows with velocities up to 50 kms"' 
are obs erved in H2O 1 1 1 -Onn, 2| 1 -- 2n2, lio-loi and 3i2-3o3 emis- 
sion ( [Emprechtinger et al. 2010 1 and up to 70 kms"' in CO 



H-6.1 kms"' H2O absorption ( [Emprechtinger et al.|2010| and the 
H-6kms"' OH component to within the 0.7 kms" velocity res- 
olution of the ATCA observations ( [Brooks & Whiteoak[[2001| l. 



7=2-1, 3-2 , 4-3 and 7-6 0achiller & Cemicharo|1990|[Leunni[ 
|et al.[[2006l l. In our HIFI observations of the CH transition, we 
detect neither blueshifted gas nor outflows to match the H2O, 
CH3OH, NH3, CO and H20^, in neither emission nor absorp- 
tion. The absence of CH in both blueshifted emission and out- 
flows indicates that CH is destroyed in shocked environments. 
We conclude that CH traces molecular material that has not yet 
been swept up by outflow activity related to star formation. 

5. Conclusions 

We use CH as a tracer of cold, quiescent, molecular gas in both 
the star-forming core of NGC 63341 and diffuse clouds along the 
line of sight. We derive a very low CH abundance of 7x10"" in 
the hot core, consistent with CH depletion in chemically active 
dense cores. Evidence is found of four absorbing components 
composed of diffuse gas. From comparison with other molecular 
absorption signatures in the literature, we conclude that clouds 
A and D are related to the NGC 6334 complex, and that clouds B 
and C are unrelated foreground clouds (Table[T]). Differences be- 
tween the four clouds in terms of counterparts in other molecular 
tracers are ascribed to local variations in density and dynamics. 
Dynamical stirring of the gas is hypothesized to be the cause of 
the absence of CH in the outflow component, which is detected 
in other molecular tracers. 
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